How far away is the next earthquake? 
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Spatial distances between subsequent earthquakes in southern California exhibit scale-free statis- 
tics, with a critical exponent 5 ~ 0.6, as well as finite size scaling. The statistics are independent of 
the threshold magnitude as long as the catalog is complete, but depend strongly on the temporal 
ordering of events, rather than the geometry of the spatial epicenter distribution. Nevertheless, the 
spatial distance and waiting time between subsequent earthquakes are uncorrelated with each other. 
These observations contradict the theory of aftershock zone scaling with main shock magnitude. 
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What do we know about earthquakes? Either a great 
deal or a meager amount, depending on the point of view 
and on the adopted definition of an earthquake. If an 
earthquake is defined to be the slip on a fault (or sev- 
eral faults) that produces the observed seismic wave field, 
then we have a good understanding of earthquakes 
In contrast to earthquake kinematics, no satisfactory un- 
derstanding exists of the physical processes in the litho- 
sphere that cause slip on faults and are, thus, respon- 
sible for the dynamics of earthquakes |2|. This dynam- 
ics entangles a vast range of space and time scales and 
manifests itself in a number of generic, empirical features 
of earthquake occurrence including spatio-temporal clus- 
tering, fault traces and epicenter locations with fractal 
statistics, as well as the Omori and Gutenberg-Richter 
(GR) laws (see Refs. 3] for a review). The Omori law 
states that the rate of aftershocks after a main shock de- 
cays hyperbolically in time |4(, while the GR law states 
that the size (measured in terms of the seismic moment 
M) distribution of earthquakes is scale- free |5(. 

The presence of vastly different scales has to be taken 
into account in order to interpret measurements correctly. 
For instance, to unambiguously measure the length of 
many natural objects, like the 'fractal' coastline of Nor- 
way, one has to specify the length of the ruler used 0, Q . 
Here we ask a similar question for distances between 
earthquakes, and find a similar result. That is, in or- 
der to determine unambiguously the statistics for how 
far away the next earthquake will be, one has to specify 
the size of the region under consideration. 

We focus on the distribution of spatial distances be- 
tween the epicenters of successive earthquakes in south- 
ern California. In the past, different possibilities have 
been proposed including power-law behavior |8| and q- 
exponential (cumulative) distributions However, 
none of the previous studies has systematically taken 
into account the physical extent of the region considered, 
nor the threshold magnitude for including events in the 
analysis. Both of these quantities have recently been re- 



vealed to be crucial for capturing robust, statistical fea- 
tures of waiting times between subsequent earthquakes 

HOE 

Here, we show that the distribution of spatial distances 
between successive earthquakes, larger than a threshold 
magnitude m, occurring within a given region of area 
L 2 exhibits: (1) power law behavior with an exponent 
S ~ 0.6, (2) finite size scaling as a function of L and 
(3) no dependence on the threshold magnitude m, as 
long as the data set is complete for that threshold. Our 
results also provide clear evidence that this behavior is 
not due to a random process bound to the geometrical 
structure of the collection of epicenters, but reflects the 
complex spatio-temporal organization of seismicity. We 
further argue that the exponent S encodes information 
about this complex dynamics, which appears unrelated 
to other known properties of seismicity. Hence S may be 
a new, independent exponent characterizing seismicity. 

To analyze the distribution of spatial distances be- 
tween successive earthquakes or "jumps" , we adopt a 
method proposed by Bak et al. and take the per- 
spective of statistical physics: Neglecting any classifica- 
tion of earthquakes as main shocks, foreshocks or after- 
shocks, analyze seismicity patterns irrespective of tec- 
tonic features and place all events on the same footing. 
Consider spatial areas and their subdivision into square 
cells of length L. For each of these cells, only events 
above a threshold magnitude m are included in the anal- 
ysis. In this way, we obtain a list of the spatial distances 
Ari = | ri+i — Ti\ between successive events with epicen- 
ters Yi and r^+i both in the same cell of linear extent L. 
Concatenating the lists of jumps obtained from each cell, 
a probability density function of the jumps P m? z,(Ar) can 
be measured 13]. Since both the threshold magnitude 
and the length scale of the cell are arbitrary, we look for 
robust or universal features of this distribution that may 
appear when these parameters are varied. 

For the SCEDC sub-catalog from southern California 
we study here, the reporting of earthquakes is assumed 
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to be homogeneous from January 1984 to December 2000 
and complete for events larger than magnitude m c = 2.4 
[l4| . Considering epicenters located within the rectangle 
(120.5°W, 115.0°W) x (32.5°7V, 36.0°7V) gives N = 23374 
events with magnitude greater than or equal to m c . The 
GR law for the cumulative distribution of earthquakes 
larger than magnitude m is P>(m) ~ iQ-bm w j iere ^ e 
seismic moment M oc 10 1,5m , and b = 0.95 for this col- 
lection of events. 

For m > ra c , we find that the jump distribution is 
described by the finite size scaling (FSS) ansatz 



P m , L (Ar) = 



f(Ar/L) 
L 



(1) 



where the scaling function f(x) decays as x with 5 ~ 
0.6 for x < 0.5 as shown in Fig.EJ For x > 0.5, it decays 
extremely rapidly since the finite cell size requires that 
f(x) = for x > y/2. These observed results have far 
reaching implications. 

First, for any L the cutoff sets in at r ~ L/2 and, 
hence, scales trivially with L. The appearance of FSS 
precludes the existence of any other length scale over the 
range where FSS holds. Thus no physical length scale 
exists in the range from 20km to ~ 500km, in contrast to 
the theory of aftershock zones . According to this the- 
ory, main shocks generate aftershocks within finite after- 
shock zones, whose extent is comparable to the rupture 
length l r = 0.02 x 10 5 m km of the main event 15]. This 
implies that the distance between subsequent aftershocks 
would be limited to the size of the largest aftershock zone, 
which is less than 90km for the catalog analyzed here. As 
in Ref. we find no evidence for this length scale. 

Further, as S is unambiguously less than one, the dis- 
tribution P m5 z,(Ar) becomes non-normalizable for large 
L. Extrapolating our results, the finite size of the earth 
may play an important role in the definition of distances 
between subsequent earthquakes. Finally, P mjj £,(Ar) does 
not significantly depend on m for m > ra c , though the 
number of included earthquakes varies considerably with 
magnitude. 

Although the distribution of jumps reflects a dy- 
namical property of seismicity, the particular form of 
P m? L(Ar) could be determined by the geometrical struc- 
ture of the collection of epicenters. A simple test can be 
made by randomly rearranging or 'shuffling' the tempo- 
ral sequence of events, while holding the magnitude of 
the events and their epicenters fixed. The distribution 
of distances between subsequent events in the shuffled 
catalog is very different from the original one, as shown 
in Fig. For the shuffled catalog, the distribution of 
jumps does not decay with an exponent 5, but it actually 
increases as a power law with exponent S s h u f — —0.14. 
This can be understood from the fact that the randomly 
rearranged ordering gives an estimate of the distribu- 
tion of distances between any two earthquakes within the 
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Figure 1: (Color online) Rescaled distribution of jumps, 
P mjJ c(Ar), for different values of m and L. L — 448.5km 
corresponds to the full area where there is just one cell. Val- 
ues of Ar < 2km have been discarded due to uncertainties 
in estimating the epicenters' locations. The solid line with 
exponent 5 — 0.6 is shown as a guide to the eye. 



same cell, which is also shown in Fig. O The critical ex- 
ponent of the cumulative distribution is, by definition, a 
measure of the correlation dimension D2 of the epicenter 
distribution [l7| • O ur findings imply a fractal dimension 
D 2 = 1 — S s huf ~ 1-14, which agrees, within statistical 
error, with the value obtained in [l^ . 

The comparison between the distribution of jumps 
and the distribution of distances between epicenter pairs 
clearly proves the dynamical origin of a non-trivial S. In 
particular, this exponent is not simply related to the cor- 
relation dimension. Yet, D2 is also independent of the 
threshold magnitude of the earthquakes considered and 
the size of the region studied 18]. The latter fact and 
the appearance of FSS in all distributions of Fig. [2J are 
inconsistent with the existence of any preferred length 
scale, other than the cell size chosen by the observer. 

The "propagation" of seismic activity is not only de- 
scribed by spatial distances Ar^ but also by the wait- 
ing time AU between successive earthquakes i and i + 1. 
Although the statistics of the waiting times has been 
studied recently [HI El E3 ? the propagation itself has 
not been analyzed. Here we consider the velocities, 
Vi = Ari/Ati, between subsequent events in each cell for 
events with magnitude above a threshold m and com- 
bine them, as before, into a probability density function 
Pm,L(v). (Note that P m ,L(v) is not directly related to the 
controversial and debated subject of aftershock diffusion, 
which refers to the expansion or migration of aftershock 
zones with time. See Ref. 19] for a review.) 

Figure |3| shows that Pm,L(v) ~ ^ _r? with 77 w 1.0 for 
intermediate v. The cutoff at large v is determined solely 
by L, due to the fixed temporal resolution of At > 60 sec 
implying v max = V2 ■ 100km -GOh -1 w 8485km/h. The 
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Figure 2: (Color online) The jump distribution, P mj z,(Ar), 
for m = 2.4 and L — 448.5km. (Orange) circles corre- 
spond to the distance between successive earthquakes Ari = 
ri+i —Ti\ as in Fig. Q] (Red) squares correspond to An = 
r cr(i+i) — r <r(i) | where a is a randomly chosen permutation of 
the integers < i < N, giving a shuffling of the catalog as 
described in the text. The (black) triangles correspond to the 
distance between any two earthquakes in the same cell. The 
solid line is a fit to the latter distribution from 2km to 200km 
with a critical exponent of 5 s h u f = —0.14 ± 0.03. 
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Figure 3: (Color online) The velocity distribution, P m ,z,(v), 
for L = 100km and different values of m: (Blue) crosses 
and boxes correspond to the original data with v% — 
An/Ati. (Red) triangles correspond to a randomization 
Vi — Ar a (i)/Ati, where a is a randomly chosen permuta- 
tion of the integers < i < N. Note that Ari > 2km and 
Ati > 60s for all i to avoid any bias due to the uncertainties 
in the estimates. Inset: Rescaled distributions, using the GR 
exponent b = 0.95 and n — 1.0. 



cutoff at small v depends on m. The inset of Fig.Elshows 
that its position scales as 10 _6m . If we ignore, for the 
moment, the variation of distances Ar, then P m ,L(v) — » 
Pm,L(l/ At). The latter distribution is obtained from 
the distribution of waiting times P m ^(At) in 0, 0, 
EH- Since P m , L (l/At) d(l/At) = P m [ L (At) dAt, in this 
approximation the exponent rj = 2 — a, where a ~ 1 is the 
exponent characterizing the distribution of waiting times 
for intermediate arguments. Furthermore, the cutoff at 
small v would be controlled by the cutoff at large t in 
the waiting time distribution where the behavior crosses 
over from a power-law regime with exponent a ~ 1 to a 
faster decay at t cutoff - 10 6 m HllJIli- 

If the statistics of the waiting times and the jumps are 
independent, then P m ,L(v) will only reflect the statistics 
of the waiting times over a range of velocities. This is due 
to the fact that the distribution of waiting times is much 
broader (approximately seven orders of magnitude) than 
the distribution of spatial distance (approximately three 
orders of magnitude) for our data. Indeed, the particular 
combination of spatial and temporal distances between 
successive earthquakes is largely random. As shown in 
Fig. 3, the estimate of P m ,L(^) does not change signif- 
icantly if Vi is given by V{ = Ar a ^\/Ati where a is a 
random permutation of the integers < i < N. This 
clearly proves that waiting times and jumps are indepen- 
dent of each other. Additionally, using the same values of 
m and L as in Fig. 03 the spatial and temporal distances 
between successive earthquakes are almost uncorrelated 



as measured by the cross correlation 

r = ((Ar, - (Ar,),) • (At, - <At i ) i )) i _ Q QJ 
y/((An - (Ar,),)2), % /((At i - (At,),)2), ~ ' ' 

where (. . . )i denotes an average over events i. 

Considering the observation that the critical exponent 
6 depends on the temporal order of events and is not 
determined solely by the geometry of the set of epicenter 
locations, together with the observation that the waiting 
times and jumps are largely uncorrelated and apparently 
independent of each other suggests that the exponent 
5 may be a new, independent exponent characterizing 
seismicity. 

The lack of correlation between waiting times and 
jumps could be interpreted as an indication that the af- 
tershock decay rate at all distances is the same, and could 
therefore have implications both for models of aftershock 
diffusion and for the rate and state/static stress friction 
model of aftershock triggering 19]. However, our analy- 
sis is not based on the distinction between main shocks or 
aftershocks. This distinction is relative 0,[2£j. In fact, 
there is no unique operational way to distinguish between 
aftershocks and main shocks and t hey are not caused 
by different relaxation mechanisms 0,[23. Besides, such 
a classification may not always be the most useful way 
to describe the dynamics of seismicity. One could never- 
theless consider the possibility that the scaling region in 
Fig. [I] can be attributed entirely to aftershocks. This 
would require that aftershock sequences dominate the 
statistics during the period and magnitude range con- 
sidered. According to the traditional classification and 
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using aftershock zone scaling with main magnitude, the 
maximum distance between aftershocks would be deter- 
mined by the largest events, namely the m = 7.3 Landers 
earthquake and the m = 7.1 Hector mine earthquake. 
Thus, this distance should be less than 90km. However, 
we find no break or change in scaling behavior for larger 
distances extending all the way up to the size of the re- 
gion consider, of the order of 500km. 

Our results are strikingly different from earlier results 
by Ito [3 and Abe and Suzuki 9] for California, who ex- 
amined similar catalogues over a similar time span but 
did not take into account the length scale of observa- 
tion. The latter authors also included earthquakes with 
magnitude as low as 0.0, and (as in 8]) found a very 
different jump distribution, implying that their analysis 
suffered from the incompleteness of the earthquake cat- 
alog at small magnitudes. Although the distribution of 
jumps is independent of the threshold magnitude as long 
as the catalog is complete, the marked difference between 
our results and previous ones shows that the complete- 
ness of earthquake catalogs is crucial for obtaining robust 
and accurate results. 

As pointed out by Corral the (long-term) mea- 

surement of the distribution of waiting times between 
subsequent earthquakes — as in Ref. [HI — involves 
an average over regions and time spans with widely dif- 
ferent rates of seismic activity. If, instead, statistics are 
measured in regions and time spans with a stationary 
rate of earthquakes, a different distribution is obtained. 
A similar situation occurs for our measurement. The 
universal law we find for spatial distances between sub- 
sequent events holds for data sets where the rate of earth- 
quakes is heterogeneous, i.e., for rather long time spans. 
Analysis at short time scales in the stationary regime, as 
described in Ref. will be explored in a future work. 

To summarize, we have shown that the distribution 
of spatial distances between successive earthquakes in 
southern California exhibits finite size scaling with a non- 
trivial power law exponent S ~ 0.6. Thus, in order to 
specify the statistics of distances between earthquakes, 
one has to define the length scale of observation, since 
no intrinsic length scale exists. This observation is in 
sharp contrast to the theory of aftershock zone scaling 
with main shock magnitude The exponent S has 

a dynamical origin, but the distances and waiting times 
between subsequent events are found to be independent 
of each other. This implies that the complex dynamics 
of seismicity has a self-similar hierarchical structure in 
space and time, consistent with the hypothesis that it is 
a self-organized critical phenomenon |24j. Our findings 
can be used as benchmark tests for models of seismicity. 
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